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Abstract 

We investigate the effects produced on the Ught scattering spectrum by the anisotropic diffusion 
of impurities (dye) in a compressible nematic solvent. This spectrum is calculated by using a fluc- 
tuating hydrodynamic description when the system is in both, a fully thermodynamic equilibrium 
state and in a non-equilibrium steady state {NESS) induced by a dye-concentration gradient. 
In the former state, the isotropic pre-transitional as well as the nematic phase of the solvent are 
considered. We find that the equilibrium spectrum is symmetric (Lorentzian) with respect to the 
frequency shifts, but anisotropic through its explicit dependence on the ratio of the parallel and nor- 
mal diffusion coefficients of the dye. The values of these coefficients were taken from experimental 
measurements of diffusion of methylred and nitrozo di-methyl aniline in a MBBA solvent. We find 
that the compressibility of the solvent increases the maximum and the width at half height of the 
Rayleigh peak, with respect to the incompressible case jlj]. This increase varies between 12% and 
25%, respectively, when the impurities concentrations is the range of 1% - 5%. The NESS induces 
a coupling between the concentration fluctuations of the dye and the hydrodynamic fluctuations 
of the solvent. In this case the compressibility effects may increase the maximum and the width 
of the central peak up to 25%, for values of the concentration gradient four orders of magnitude 
smaller than those considered in the incompressible case. This result indicates that compressibility 
and mode coupling effects may be significant and that they might be detected experimentally. On 
the other hand, for the nonequilibrium Brillouin part of the spectrum we find that the intensities of 
the sound propagation modes are unequal and one of the peaks shrinks in the same amount as the 
other increases. This asymmetry increases linearly with the magnitude of the solute concentration 
gradient. The maximum difference between the nonequilibrium and equilibrium contributions to 
the Brillouin spectrum for various values of the external gradient is also estimated. However, in 
all cases we find that the amplitude of the Brillouin peaks is several orders of magnitude smaller 
than the central peak. Therefore, although nonequilibrium effects do produce an asymmetry, our 
theoretical analysis indicates that these effects are too small to be observed experimentally. 

PACS numbers: 77.84.N, 61.30G, 42.65 
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I. INTRODUCTION 

In a previous paper hereafter referred to as / we have analyzed the effects pro- 
duced by the anisotropic diffusion of impurities in an incompressible nematic solvent. By 
using a fluctuating hydrodynamic approach we calculated the central Rayleigh peak of the 
spectrum when the solvent is in both, a fully thermodynamic equilibrium state and in a non- 
equilibrium steady state (NESS) induced by a dye-concentration gradient. The nonequilib- 
rium states considered in I were close to equilibrium and the driving concentration gradient 
was taken into account through a local version of the fluctuation-dissipation theorem for the 
stochastic current of the impurities. In the present work we extend the analysis of I in two 
aspects. First, instead of introducing spatial inhomogeneities in the fluctuation-dissipation 
theorem, we consider mode coupling terms relating the concentration fluctuations of the 
solute and the orientation and velocity fluctuations of the solvent. Secondly, apart from 
the central peak, in the present work we also calculate the Brillouin part of the dynamic 
structure factor of the impurities and analyze the mode coupling contribution to the full 
spectrum of the suspension. For this purpose the compressibility of the solvent is considered 
explicitly in its equations of motion. 

We find that in equilibrium the compressibility of the solvent increases the maximum 
and the width at half height of the Rayleigh peak, with respect to the incompressible case, 
by amounts that vary between 12% and 25%, respectively, as a function of the impurities 
concentrations in the range of 1% - 5%. In the NESS induced by the concentration gradient, 
these features of the central peak may increase up to 25%, even for values of the concentration 
gradient which are four orders of magnitude smaller than those considered in /. This result 
indicates that compressibility and mode coupling effects may be significant and that they 
might be detectable experimentally. 

On the other hand, for the Brillouin part of the spectrum in the NESS we find that the 
intensities of the sound propagation modes are unequal and one of the peaks shrinks in the 
same amount as the other increases, a behavior that is also predicted and observed for a 
simple fiuid [sl, We find that this asymmetry increases linearly with the magnitude 
of the solute concentration gradient. The maximum difference between the nonequilibrium 
and equilibrium contributions to the Brillouin spectrum for different values of the external 
gradient is also estimated. However, in all cases we find that the amplitude of the Brillouin 
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FIG. 1: Schematic representation of a plane homeotropic cell with a constant thermal gradient 
along z direction. The inset shows the scattering geometry.The scattering angle is d. 

peaks is several orders of magnitude smaller than that of the central peak. Therefore, al- 
though nonequilibrium effects are present in the above mentioned asymmetry, our theoretical 
analysis indicates that these effects are too small and difficult to observe experimentally. 

II. MODEL AND BASIC EQUATIONS 

As in /, we consider a dilute suspension of noninteracting impurities diffusing through 
a thermotropic nematic liquid crystal solvent, as depicted in Fig. 1. If the suspension is 
sufficiently diluted, the dynamics of the impurities does not disturb appreciably the state of 
the nematic and it may be considered to be in an equilibrium state defined by a temperature 
To, a pressure Po, a vanishing velocity field ~v° = and a uniform director's orientation 
n° = (0,0, 1), corresponding to the homeotropic configuration shown in Fig. 1. X 

We shall only consider nonequilibrium states corresponding to a stationary concentration 
field of impurities defined by 

c" (7*) = Co + ~r^ -li, (1) 

where Co is the mean concentration of impurities and li = Vc is the uniform concentration 
gradient in the y — z plane and whose direction is specified by ip. For future use it will be 
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convenient to recast ([T]) in the more convenient form 

w— >\ 1- sin("g'-"r') 
c' (r) = Co + a hm — -, (2) 

q — ►o q 

where a = \~a\, ~q = ^Vc is an auxihary vector of variable magnitude and parallel to ~a . 

If no chemical reactions occur between the impurities, their total number is conserved 
and their local concentration density, c(T^,t), obeys the continuity equation 

I + V.J. = 0, (3) 

where Jj (V, t) is the flux of the diffusing particles, which for an uniaxial nematic is of the 
form 

ji (r, t) = -Dij (r , t) v,c (r, t) + c (t^, t) Vi (r, t) . (4) 

The first term on the right hand side is the usual Pick's law contribution where Dij (T^, t) is 
the diffusion tensor of the suspended impurities; the second term represents the convective 
diffusion of the impurities, where Vi {'r^,t) denotes the velocity field of the solvent. For an 
uniaxial nematic Dij has the standard form 

Dij (r, t) = D^5^j + {D\\ - D^) m (r, t) n, , t) , (5) 

where D\\ is the diffusion coefficient parallel to the director n{'r^,t) and D_l is the corre- 
sponding coefficient in the perpendicular direction. Da = -Dy — -D_l is the corresponding 
diffusion anisotropy. Usually, the diffusion of small particles dissolved in a nematic solvent 
is such that the diffusion parallel to the director is faster than perpendicular to it; as a 
consequence, the ratio D\\l Dx_ seems to be independent of the actual shape of the diffusing 
molecules jst, [gI, [3]. Using Eq. ([5]), the diffusion equation for c(^,t) turns out to be 

d 

— c(r,t) = Diy^c + DaTiirijViVjC 

+Da {niVifij + njVirii) Vjc 

-ViViC + cViVi. (6) 

Note that the compressibility of the solvent is taken into account by the term ViVi ^ 0, 
which was absent in /. 

If a fluctuating mass diffusion current, Jf (l^,t), is introduced into this equation, the 
concentration fluctuations 6c{l^,t) = c(T^,t) — (^), obey the linearized equation 

^<5c = (DxVi + D|| V^) 6c + - c'VM 

+Da {a^VJn, + a.Vidrii) - V^Jf. (7) 



Jf {'r^,t) is a Markovian, Gaussian, stochastic processes with zero mean (^Jf {l^,t)^ = 
and whose correlation is assumed to obey a local equilibrium version of the usual fluctuation- 
dissipation relation j^, 



If we define the Fourier transform of an arbitrary field / ("r^, t) by 
Eq. ([7]) reads 

6c (^k ,u?j = G (^k ,u?j [—iCoki6vi + ~a ■ V-^ {ki5vi) — ai6vi 
+iDa {kzaiShi + azkiShi) — ikiJ[ 



where the propagator G [ k ,Ljj is given by 

G (t, u?j = {-iuj + D^kl + D\\kl) 



(8) 



(9) 



(10) 



(11) 



with k"^ = kl + ky and where V-^ denotes the differential operator in k -space, V-^ = 
{d/dkx,d/dky,d/dkz). Similarly, Eq. ([H]) becomes 



jF(J,uj]jf( k\J 



2 {2-Kf DijS {uj + uj')\co5{k + k' 



a 

+ 2^ 



5ik + k' -5ik + k' + ^ 



(12) 



A. Impurities Structure Factor in Equilibrium 

The impurities dynamic structure factor in equilibrium, S'^" ( , ) , is obtained by setting 
^ = in Eqs. ^ to yield 

Sc (^k ,u?j Sc k , —CO 

sfCk,u) + s:^(~k,uj 



(13) 



S"' k,u 



2(27r 
- G{ k,u 



Gi k,uj 
2 



5^ (0) Cohk^Di 



c^kikj (6vi[ k ,uj] 6vj { — k , —uo 



The first term is the contribution due to the stochastic current J[ and the second one 
is the contribuiton arising from the dynamics of the nematic solvent through its velocity 



correlation functions. Note that in contrast to /, the equihbrium correlation function of the 
solvent velocity fluctuations appears explicitly in this equilibrium property of the solute. 
This correlation will now be calculated from the fluctuating hydrodynamic equations for the 
solvent. 

The hydrodynamic state of the nematic is described in terms of the pressure p{~r^,t), 
temperature, T{~r^t) and the velocity ~v{~r^t) fields, and the unit vector defining the local 
symmetry axis (director), n(l^,t). To calculate (^5vi{^k 5vj k ,—u?j^ and other 
correlation functions of the solvent that will be appear below, it is convenient to separate 
the state variables into two independent sets, namely, transverse and longitudinal variables 
with respect to the plane h° — k , j^, [lO]. The former set is {Sfii, 6vi} with 



Shi = kj_^h° 



5vi = kl^h'' 



k X 6 n 
k X 6~u 



while the longitudinal set is < 611,3, 6v2, ST, Sv^, 6p > with 



5v'^, = k kj_ k X 



k X h° 



5lf 



5v3 = k ^ k ■ 6~n 



and 



(14) 
(15) 

(16) 
(17) 

(18) 



= k ^ k ■ 5~ri , 

Substitution of these definitions in Eq. (fT3l) . evaluation of the resulting expression at k — 
— k , uo' = —uj and use of the explicit expression of the propagator G {^k leads to the 
following equilibrium structure factor 

ijJoik) ( Cok^ 



S"'{k,iu) = 2(27r)'5^(0)co 



1 



ujd{ k 



6V3 [ k ,uj] 5f 3 [ — k , —uo 



with 



uJoik] =D^kl + Dnkl 



(19) 
(20) 



B. Nonequilibrium Impurities Structure Factor 



Let us now consider the effect produced by the concentration gradient Vc. The nonequi- 
librium part oi S (^k , cu^ contains contributions arising from jf, which is not coupled with 



Svi, Srii, and three different contributions arising from the dynamics of the nematic solvent 
which are expressed as director, velocity and cross director- velocity correlation functions, 
that is, 

(t, 00) = ST' (t, + SZ" (t, 00) + SZ' , + SZ' (t, . (21) 



Sjyk,uj is the nonequilibrium contribution arising from J/* due to the assumption of 
the validity of the local version of the fluctuation dissipation theorem (jH]), which is given by 
Eq. (30) in /, 



(5c I k ,u] Sc \ — k , —ijj 



neq 



2 ( 6c \ k ,uj] Sc \ — k , —oj 



eq 



IVcl ~ \ 2 

X a; — — {D\\kz cos Tp + Dx^k^sinilj) G yk ^ojj 

Furthermore, since longitudinal and transverse fluctuations are uncoupled, 

2 



-Di 



G[k,u 



k1a\ [ k ] ( 5ni [ k ,uj] 6ni { — k , —cu 



+ 



2k,a3 ( k ] + {k^- ] a2( k 



X ( 6n3 [ k ,u) 6n3 { — k , —u 1,1, 



(22) 



(23) 



a k,u 



Da 

+2cok 

X Re < ( (5n3 ( k ,uj) 5v3 ( — k , —uj 



—2kza^ ( k j Im | ^ 5ni ( A; , ) 5vi 

kz 



k , —uj 



2kza3 ( ) + I 1 - ( ^ ) I k^a2 ( k 



+2a2 k 



2kza3 ( k 



1 _ ( _± ) I k^a2 ( k 



X 



3^ |^(^'y2 (^k ,u?j 



5n-^ i — k , —UJ 



(24) 



SZ'{k,uj 



G k,Lj 



0.1 [ k ] ( Svi [ k ,uj] Svi [ — k , —to 



+ a\ [ k ] \ UV2 \ ,LU j UV2 



5v-),[ k ,uj] 5v2 { - k ,-uj ^ ^ 
-2coka2 ( k ) Im <j ( 5v2 ( k ,u) 5v3 [ - k , -cu) ^ 1 1 



(25) 
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where a^\^k j, jj, = 1,2,3, represent the transverse and longitudinal components of the 
concentration gradient, defined in a similar fashion as 6h^ and (Eqs. f|T4l) - f|T8l) ). It should 
be pointed out that Eqs. (!23|) - (!25l) show that in contrast to equilibrium, in NESS the density 
gradient introduces a coupling between the concentration fiuctuations of the solute and the 
velocity and orientation equilibrium fiuctuations of the solvent. These contributions should 
be calculated by first evaluating the required correlation functions of the solvent. 



III. SOLVENT EQUILIBRIUM CORRELATION FUNCTIONS 



Let us recall that the hydrodynamic fiuctuations of a thermotropic nematic evolve on three 
widely separated time-scales corresponding to the the relaxation of orientational, visco-heat 
and sound modes, respectively, 13], [ll]. These relaxation times are such that Torientation ~ 
ulKk"^, Tyisco-heat ~ pCp/ nk'^ ~ p/uk"^ and Tsound ~ ^ / Cgk , where u denotes any of the 
nematic's viscosities, K is the elastic constant, Cp denotes the specific heat at constant 
pressure, k is the magnitude of any of the components of the thermal conductivity tensor 
and Cs is the isentropic sound speed of the nematic. For values of k corresponding to a 
hydrodynamic description and for typical values of the material parameters of a thermotropic 
nematic [12|], the following relation holds 

'^orientation "^visco-heat > Tsound- (26) 

By estimating the order of magnitude of the elements of the hydrodynamic matrices of the 
time evolution equations for the nematic's fiuctuations, which are given in Ref. [l7] by Eqs. 
(24), (25), (38)-(42), it is possible to identify the following groups of variables {6hi,6h3}, 
l^-fii, 5-02, 5r|, {5v3,Sp} as slow, semi-slow and fast, respectively. The wide separation 
between these time-scales will now be exploited to eliminate the faster variables from the 
general dynamical equations obtaining, thus, a reduced description in which only the slower 



variables are invo 



ved. For this purpose we use the time-scaling perturbation method devel- 



oped in Refs. [13|], which allows to find a contracted description in terms of the slow 
variables only. The corresponding reduced dynamical matrix will be constructed by a per- 
turbation procedure, where the perturbation parameters are the ratios Tmsco-heat / ^orientation 
and Tso«rid/'i^«sco-/ieat,[l^- Usiug this formalism it can be shown that in the slow time-scale. 
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the director fluctuations Srii and Sris obey the stochastic equations 



iu5h^ = -Unfj, {^k^ 5n^ - (Jn^, /x = 1, 3, 



where the fluctuating terms cr„^ are 
1 



0"ni = 



2 V2ki + y^kl 



;i + A) + (1 - A) kl 



2 '"y'^j^yj 



2 z/gA;^ + 2 (z/1 + z/2 - z/g) A;2 + u^k^ 



k kj ^-j^j k^, k-i kj 



(27) 



(28) 



(29) 



Here Tj and Sjj denote the stochastic components of the director's quasi-current and the 



stress tensor, Uijki, which obey the the fluctuation-dissipation relations 13] 

{T^ij (7", t) Sh (7", t')) = 2kBToUijkiS (7" -1^)6 (t' - t) , 



(T,(r,t)T,(r',t')) 



7i 



(r' - r) <5 (f - 1) 



(30) 
(31) 



where fc^ is Boltzmann's constant, 71 is the orientational viscosity coefficient. The viscosity 
tensor is 



M^jiSik+5ii5jk) + 2(z/i z/2 - 2z/3)n>°nfcn° 

+ (z/3 - z/2)(nX'^ifc + + n°n°Jji + n^X^jk) 



(32) 



where z^i, 1^2 and z/3 are three shear viscosity coefficients and z^s and z/4 — z/2 denote two bulk 
viscosity coefficients. The quantity 



= Sij - n°n° 



(33) 



is a projection operator and 5ij denotes a Kronecker delta. In the above equations we have 
also used the following abbreviations 

l7i(l + A)'A;2" 



I k 



1 

7i 



{K2ki + K^e;) 



1 + 



UJn-A I k 



4 z/2/i;i + z/3A;2 
7i[(l + A)fc2 + (l-A)A:i]^ 



1 {K,kl + K,kl) 1 + 1-^^^, + 2 (z/, + z/2 - z/3) + u,kt 



(34) 



(35) 



where i^i, K2 and ii's are, respectively, the splay, twist and bend elastic constants and A is 
a non-dissipative coefficient associated with the director's relaxation. 
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From Eq. (31) in Ref. 17|, we obtain the reduced equation for the semi-slow transverse 
variable, 6vi, 

— iu6vi = —ujy-^ ^ A; j 6vi — 5"^^ , (36) 
where the stochastic force term a^^ is given by 



a. 



Pok. 



(37) 



and the reduced characteristic frequency u^-^ ( A; ) is 

1 



Po 



(38) 



Similarly, from Eqs. (39)-(42) in Ref. 17|, we obtain the corresponding equations for the 
semi-slow longitudinal variables 5v2 and (5T, 

^v2(yk^ 
ujt [ k 



— lUJ 



6V2 

Sf 



)l 


( 6V2] 


l-l 


' 5-^2 1 




iSf) 







(39) 



with 



ft/,, Zj , 



kikj^ij 

k\k 



ax = i kiQi. 

Po^^p 



(40) 
(41) 



Here Qi is the stochastic heat flux which satisfles the fluctuation-dissipation theorem 

{Q^ (^, t) Q, (r', t')) = 2kBT^K,,6 (r' - r) 6 {f - t) , (42) 
Kan°n° with the heat conductivity tensor and uj^^ik], 



uj^ [ k ] are given by 



uJv2 I k 



1 



h(A:,' + fcl)+2(z/i + Z/2-^^3)^'fci] 



(43) 



Pok"^ 

uj^ (^"fc ) = D^fci + Dlkl, (44) 

where = i^j^/ PoCp-, Dj^ = K\\j p^Cp^ stand for the thermal diffusivity coefficients of the 
nematic along the directions parallel and perpendicular to h° . 

Analogously, from Eqs. (60) in Ref. [l7|, we obtain the reduced equation for the fast 



variables 5f 3 and 5p 



6p 



( 



^vz ( k 



ic^k I i— 
IP 



1/2 



1/2 \ 



\^csk[^^) {^-l)uT[k) j 
11 



5p 



(45) 



which is vahd in the fastest time scale. The stochastic noises a^^ and dp are defined by 



(7^,3 — -kikjTiij, 



(7p — ICg 



(7-1) 



1/2 



(46) 
(47) 



where we have used the definition 



+ (2z/i + + 2z/5) fc^] . 



(48) 



It is essential to stress that the dynamic equation Eq. (!27j) is correct in the slowest time 
scale, that is, for times of the order of Torientation- Similarly, Eqs. ( l36l) and ( l39i) are valid for 
times of the order of Tmsco-heat and Eq. (H5l) describes the dynamics of the fast variables in 
the fast time-scale characterized for times of the order of Tgound- As a first approximation we 
extrapolate Eqs. fl36|) . fl39|) and fH5l) to the slow time-scale, in order to calculate the required 
nematic's correlation functions. Solving Eqs. (127|) . (l36ll . (l39ll and (145|) for (5n^ and 6v^, using 
the definitions of the stochastic terms dn^i, o'^^, 5"t, 5"p and using the fiuctuation dissipation 
relations (I30l), (EI]), (02]), we arrive at 



5ni ( k ,u) 6ni { — k , —lu 



0L\ ( k 



71 + ( A; 



(Jn-? ( , cij ) 5ns ( — , — cij 



^2 as A; 



(5f 1 I A; , cij ) 5t>i I — A; , — 



bv') \ k bv-y \ — k , —uj 



71 uj^ + ul (t) 
^ A; ^ 



Pou^ + Ul { k 



g idv^i k 



Po UJ^ + ul i k 



(49) 



(50) 



(51) 



(52) 



and 



5{).s { k ,uj] 6v^ \ — k , —UJ 



Po 

X < ti;„ ( k 



UJ + Cskf + r2 ( A; ) {uj- Cskf + V^[k 



+ {-i-lful ( k 



+c;k^ (7 - l)cjj, k 



(53) 
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where e = 2(27r)^ (5^(0) keTo 



7i(l + A2)^' 



ai( k 



1 + 



asik) = 1 

r ( k 



A{u2ki + p^kiy 

71 [(1 + m + (1 - \)kif 

4 [z/3 (fcl + ki) + 2 (z/i + z/2 - ^^s) A;iA;2] ' 
1 



+(7-1)^ f k 



(54) 

(55) 
(56) 



r 1^ A; j is the anisotropic sound attenuation coefficient of the nematic. To arrive at the 
previous correlation functions we took into account that Tj, Sj^ and Qi are not correlated, 
and that for typical values of the material parameters of a thermotropic nematic the relation 
Cgk ^ ujy^ {^k^ ,uot (yk^ ^ uonfi ^ ) > which are equivalent to ( l26l) . holds. 
Following the same procedure described above, it can be shown that 

,2 



Im < ( (5fii ( k ^uo] 6y 



k , —uj 



Po ul, ( k ) 



X 



u^ + ulAk uj^ + ulAk 



(57) 



Re < ( (5n3 [ k ,uj] I — A; , —oo 



Im 



{(*v.(- 



k ,uj] 5nq \ — k , —uj 



e k 



[{l + X)kl + {l-\)kl] 



2po k^ 

(7 - 1) u;r ( k 



^ 1 

cj^ + ivlJk^ 2 

r (li 



r A; 



(a; + csky + r^k 



(uj - c,kf + T^(k 



(58) 



e (1 + A) A;,^ + (1 - A) fcj fcV 



^^2 ( k 



X 



c.2 + ^2^( ) u^ + ul ( A; 



(59) 
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(60) 



These expressions determine the required nematic's correlation functions, Eqs. fl23l) - fl25l) . 
IV. RESULTS 

In / we showed that in equihbrium, the main contribution to the dynamic structure 
factor of the impurities is a central Rayleigh lorentzian peak. However, when compressibility 
effects are considered, the equilibrium dynamic structure factor of the impurities involves 
the equilibrium auto-correlation function of the fluctuating component of the velocity along 
k . In previous work we have shown that this correlation function contains information 
about the propagating sound modes of the nematic which gives raise to its Brillouin peaks 
[isl . Therefore, it can be expected that the light scattering spectrum of the impurities will 
also show these features and it should also exhibit two Brillouin-like peaks. For this reason, 
hereafter we will only consider the behavior of the dynamic structure factor, S (^k ,u?j, for 
frequencies close to = and u = Cgk. 



The evaluation of the different contributions ofSyk,ujj can be simplified by considering 
the order of magnitude of the involved material parameters. The diffusion coefficients of dyes 
in a thermotropic nematic are of the order of D ~ 10^^ cm^ j^, while for a typical room 
temperature thermotropic we have po ^ 1 g cm~^, Cg ~ 10^ cm v ~ 10^^ poise, 
~ 10^^ cm? and K ~ 10^^ dyn. At low concentrations (5%) Co ~ 10^° cm~^. 
This implies that the characteristic diffusion time of the impurities is much slower than 
the corresponding one to the director relaxation and, therefore, to all the other dynamic 
processes, ^ A; j ^ ujd ( ^ )• Therefore, by inserting Eqs. (I49l) - (l60l) into Eqs. (123!) - (!25|) 
and retaining only the leading terms corresponding with the previous orders of magnitude 
of the material parameters ata;~u;z)(A;) (Rayleigh peak, R) and at a; ~ c^/c (Brillouin 
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peaks B), we obtain explicit expressions for the different contributions of S yk ,u!j which 
are given below. 



A. Equilibrium Light Scattering Spectrum 

1. Central Peak 



In order to compare the relative effect of the compressible character of the solvent and 
the external concentration gradient on the spectrum of the impurities, we define the dimen- 
sionless structure factor S (^k ,u?j by 



S i k ,uj 



S(~k,uj) = (61) 



where S^^ (^k ,u?j represents the structure factor of an incompressible nematic in the equi- 
librium state. When the incompressibility condition is implemented in Eq. (fT9l) we obtain 



S2 {k,uA^2 {2ixf 5' (0) c„ (62) 



From Eqs. (|T9|) and fl53l) it follows that the equilibrium dynamic structure factor of the 
impurities for small frequency shifts, i.e. u;~u;£)^A;j,isa Lorentzian given by 

ZoksTo - 1) ( k 



(63) 



/^o'^s UJd (^k 

where the normalized frequency cOo = ^I^d (^^) been used. The second term on the 
r. h. s. represents the contribution due to the nematic's compressibility. This can be seen 
more clearly if this term is rewritten in terms of the isentropic compressibility defined by 
the thermodynamic relation Xs = ^/Poci- The value of this term can be estimated by taking 
typical values of the involved parameters. Indeed, in the case of the diffusion of two different 
dyes (methylred and nitrozo di-methyl aniline) at the room temperature in the thermotropic 
nematic MBBA at low concentrations (5%), we have 



^2 



10-\ (64) 
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FIG. 2: Normalized central peak in equilibrium, 5^ , as function of the normalized frequency 
iOo- ( — ) corresponds to a compressible nematic solvent and (- - -) denotes the incompressible 
contribution obtained in /. S'j^ is calculated from Eq. (|63p for typical values of the material 
parameters and the scattering geometry shown in Fig. 1 with ki = 1 x 10^ cm^^ and ^ = 1°. 

which implies that the compressibility contribution to the central peak may be significant 
and of the order of ~ 10%. To illustrate this effect quantitatively we consider a fixed k 
corresponding to an incident wave with ki = 10^ cm~^ for a scattering angle 6 = 1°, in the 
scattering geometry of Fig.l. If we plot S"^ for both, the incompressible and compressible 
cases as functions of Uo, we get the curves shown in Fig. 2. Note that the dynamic structure 
factor in the compressible case (continuos line) is higher and wider than in the incompressible 
situation (dashed line). For Cq = 2 x 10^^ cm~^ (diluted suspension at 1%) the relative 
differences of the height and half width at half height are 5% and 2.5%, respectively. For 
Co = 1 X 10^° cm~'^ (diluted suspension at 5%) this changes are 25% and 11% , which may 
be significant. 

2. Brillouin Peaks 

The Brillouin peaks are located at the frequencies u; ^ ± Cs/c and Eqs. (fT9ll and ( |53ll yield 
the following expression for the normalized Brillouin spectrum of the impurities in terms of 
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+ 



1 



UJn — 



+ r2 



(65) 



with To = T(k ) /ujd ( k 
Brillouin peaks is 



First, notice that the ratio of the maxima of the central and 

kBToCo 



c 



PoClTo 



(66) 



Thus, if we consider the order of magnitude of the involved parameters for the diffusion of 
dyes in a typical thermotropic nematic as above, it follows that Q ~ 10^^'^. This shows that 
the Brilloiun component of the light scattering spectrum of the impurities is negligible when 
compared in front of its central part. 



B. Nonequilibrium Light Scattering Spectrum 

We now consider the effect of the concentration gradient on the dynamic structure factor 
of the impurities when the solvent is in its nematic phase. It is convenient to introduce the 
normalized concentration gradient components by 



kiCo 



(67) 



1. Central Peak 



If we keep only the dominant terms Eqs. fl23l) - fl25l) . for in the range 10~^ < < 1, we 
obtain that at low frequencies, ^ c<j£) ^ A; j , S'!^^'^ ^ /c , j ^ S^^^ i^k ^uj^ ^ S'"^'' {^k ,uo^ 
and the leading nonequilibrium contribution to the dynamic structure factor can then be 
written in the form 



CokBTo J 2 

PoUJd k 



^2 
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uj^ii k ] 10.^2 k 



(68) 
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Thus, in the nonequihbrium state the central peak of the dynamic structure factor, Sr — 
Sr + Sl^'^, takes the form 



Sr {uJo) 



1 



1 + 



c„A;s7;(7-l)^T(t) 



ud ( k 



Cok'sT^ 2 



7^2 



^2 



UJ. 



(69) 



This result indicates that the effect of the concentration gradient increases both, the height 
and the half-width at half-height of the spectrum. The relative magnitude of the nonequi- 
hbrium contribution is measured by the function 



Sr-si 



CokBTo ,2 
Po^D 



7^2 



^2 



+ 



IjJ. 



CO, 



V2 



If we take ~ 10^" cm-^ To ~ 300 K, D ^ 10"^ cm^ s^^ z/ ~ 10-^ poise, Po 
ki ~ 10^ cm~^ and small scattering angles, ^ ~ 1°, a typical value is 



CoksTo J 2-2 inl2-2 

-k^a ~ 10 a . 



(70) 

1 gcm~^, 

(71) 



PoDuk^ 

This suggests that the nonequihbrium contribution could be significant, ~ 10%, for nor- 
malized concentration gradients as small as a ~ 10~^, which are four orders of magnitude 
smaller than those gradients used in /. In Fig. 3 we compare Sr {lUo) with its equilibrium 
component (lOo) for the diffusion of dyes in the thermotropic nematic MBBA at low 
concentrations (5%), for the following specific values of the normalized concentration gra- 
dient, Qx — tty — 0, ttz — 1 x 10"^, and for the scattering process shown in Fig. 1 with 
ki = 10^ cm~^, 6 = 1°. We notice that the spectrum becomes higher and wider when the 
concentration gradient is present (continuos line) than in the equilibrium case (dashed line). 
For the considered values of the involved quantities the increment in the height is about 
25% and the change in the half-width at half-height is ~ 12%. 



2. Brillouin Peaks 

Following the same procedure, we find that the leading contribution to the nonequihbrium 
part of the dynamic structure factor at a; ~ ±CsA; arises from 5^^^ ( k ,uj\. More specifically. 
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FIG. 3; Normalized central peak Sfi (wo) &s given by Eq. ([69]) for fci = 1 x 10^ cm'^ and 9 = 1°. (- - 
-) represents the equilibrium contribution. ( — ) denotes the dynamic structure factor in the NESS 



induced by a normalized concentration gradient Ht = {ax, ay, a z) with 

(V^ = 0°). 



0, = 1 X 10^6 



from the cross correlation ^5^3 (^k ,uo^ 5v^ ^— k , —uj^ ^ and the normalized form given by 
(l6T|l we get 



ST M 



2k,a^ U + 1 



k,A{T]{-i -l)ujT[k I ^ 




k±a2 ( A; 



^n3 ( ^ 



(72) 



where 



A 



(T) 



;i + A)A;2 + (l-A)fc^ 



z/3 {ki + ki) + 2 (z/i + z/2 - i/s) klkl 

+ (2z/i + z/2 - 2z/3 - Z/4 + Z/5) fc^] . 



[(2l/3 + z/5 - z/2 - Z/4) fcj 



(73) 



Notice that the nonequilibrium term is an odd function of the frequency and, as a conse- 
quence, the external concentration gradient induces an asymmetry in the spectrum in such 
a way that one of its Brillouin peaks increases while the other decreases in the same amount 
with respect to their equilibrium counterparts. This effect is linear in the concentration 
gradient magnitude. Furthermore, since it can be readily shown that the function does 
not varies considerably over the frequency intervals of the order of F ( A; ) around ±Csk, we 
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can make the approximation 
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(74) 



n3 



where the upper sign corresponds to the peak located at —Cgk and the lower sing to the 
peak at Csk. The relative magnitude of this effect is given by the quantity 
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Sf 



qeg 



qeq 
^B 
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2k,as + 1 




id, 



(75) 



whose significance can be estimated by taking into account the order of magnitude of the 
involved parameters and introducing the normalized gradient components according to ()67p . 
In this way we find that 

X ~ lO'^^kia cm. (76) 

This implies that the nonequilibrium contribution to the Brillouin part of the spectrum 
could be significant only for normalized gradients of the order of a ~ 10~^, which are much 
larger than those considered for the central peak. Moreover, if the angular dependence of x 
is taken into account, it turns out that this contribution is actually one order of magnitude 
smaller. 

In order to complete our analysis we now compare the normalized Brillouin component 
of the dynamic structure factor of the impurities, Sb (^k ,iuj , with respect to S"^ (^k ,u?j, 
using the same values of the material parameters as before. For instance, for a normalized 
concentration gradient with components a^,. = = 1 x 10^^, = 0, and the scattering 
process shown in Fig. 1 with ki = 10^ cm^^ and 6 = 90°, the height of the Brillouin peak 
located at u = —Cgk increases ~ 7% while the one located at u; = Csk decreases by the 
same amount, as shown in Fig. 4. Finally, we stress that both, the Brillouin component 
of the spectrum and the nonequilibrium effect on the Brillouin peaks are several orders of 
magnitude smaller than the central component and the nonequilibrium effect on this peak, 
respectively. Thus, the possible experimental observation of the effects discussed in this 
work seems to be more feasible for the central component of S ( k ,00). 
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FIG. 4: Normalized Brillouin spectrum Sr {lOq) as given by Eqs. (j65|) and ()72p for ki = 1 x 10^ 

cm~^ and 9 = 90°. ( ) represents the equilibrium part of the spectrum while ( — ) is the dynamic 

structure factor in the NESS induced by a normalized concentartion gradient with components 
cix = CLy = 1 X 10^^, = 0. 

V. CONCLUDING REMARKS 

Summarizing, by using a fluctuating hydrodynamic approach we have investigated theo- 
retically the influence of the effects produced by a uniform impurities concentration gradient 
on the light scattering spectrum of a suspension in a compressible nematic solvent. We com- 
pared both cases, when the solvent is in a fully thermodynamic equilibrium state and in a 
non-equilibrium steady state induced by a dye-concentration gradient. In the former state, 
the spectrum is symmetric (Lorentzian) with respect to the frequency shifts, but anisotropic 
through its explicit dependence on the diffusion coefficients of the dye, parallel and nor- 
mal to the mean molecular axis of the nematic. The values of these coefficients were taken 
from experimental measurements of diffusion of methylred and nitrozo di-methyl aniline in 
a MBBA solvent. Our results showed that the compressibility increases the height and the 
width at mid-height with respect to the incompressible case in amounts which vary up to 
25% for a dye diluted suspension at 5% in MBBA. 

As was discussed above, the nonequilibrium correction turns out to be several orders 
of manitude larger for the central peak of the spectrum than for the Brillouin part. The 
Rayleigh peak becomes higher and wider when the concentration gradient is present with 
respect to the equilibrium case. For the considered values of the involved quantities, the 
increment in the height is about 25% and the change in the half-width at half-height is 
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12%, as indicated in Fig. 3. The size of this effect depends on the square of the gradient 
components. 

To our knowledge, the physical situation dealt with here has not been considered in the 
literature and our model calculations yield new results that might be observable; however, 
this remains to be assessed. 
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